JOURNAL OF RESEARCH of the National Bureau of Standards— B. Mathematics and Mathematical Physics 

Vol. 65B, No. 2, April-June 1961 

Truncations in the Method of Intermediate 
Problems for Lower Bounds to Eigenvalues 

Norman W. Bazley and David W. Fox l 

(February 10, 1961) 

Two new procedures are developed for determining lower bounds to the eigenvalues of 
linear operators. The methods are based on the theory of semibounded self-adjoint oper- 
ators in separable Ililbert space. Computation of the lower bounds is reduced to the solution 
of matrix problems. The procedures have immediate application in the estimation of 
eigenvalues and eigenvectors of differential operators appearing in quantum mechanics. 



1. Introduction 

In this paper we give two now procedures for 
determining lower bounds to eigenvalues of semi- 
bounded, self-adjoint operators. These procedures, 
while phrased in the language of Ililbert space, 
result in useful computational methods for the esti- 
mation of eigenvalues and eigenvectors of differential 
operators occurring in classical and quantum me- 
chanics. Our procedures are based on the method 
of intermediate problems. 

Intermediate problems, us introduced by A. Wein- 
stein [7], 2 give lower bounds by changing the bound- 
ary conditions of differential operators. By relaxing 
the boundary conditions he obtains a solvable base 
problem that gives rough lower bounds to the eigen- 
values. He then links the base 4 problem to the given 
problem by a sequence of intermediate problems 
which can be solved in terms of the base problem 
and which improve the lower bounds. 

Years ago H. Weyl [8] pointed out that rough 
lower bounds can also be obtained by changing the 
operator. This fact was used by N. Aronszajn [11 
to construct intermediate problems for those cases 
in which a base problem can be obtained by neglect- 
ing a positive term in the operator. In the same 
paper Aronszajn systematized the method of inter- 
mediate problems for completely continuous oper- 
ators in Hilbert space. As Aronszajn points out, 
the solution of the intermediate problems in his 
formulation requires the determination of the poles 
and zeros of a meromorphic function usually given 
as an infinite sum. 

Recently H. F. Weinberger showed that it is 
possible to overcome the principal difficulties that 
arise in the intermediate problems of Weinsteint, 
Weinberger's result, which is a special case of his 
more general theory [6], gives lower hounds to the 
eigenvalues by the computation of the zeros of a 
polynomial. 

In this paper we first sketch Aronszajn's method 
of intermediate problems for changes in the operator. 
Our presentation allows for the presence of continu- 
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ous spectra and is in a form suitable for application 
to many quantum mechanical problems. In the 
second section we show how it is possible to modify 
the method so that lower bounds can be obtained by 
elementary calculations. These results give formu- 
las parallel to some of those of Weinberger [6]. A 
feature of our procedure is that it can be used to 
calculate lower bounds to the energy levels of most 
atomic systems. The last section extends another 
method of one of the authors [2]. 

Part of the work presented here was started while 
the authors were in the Institute for Fluid Dynamics 
and Applied Mathematics, University of Maryland, 
and was summarized in [3]. 

2. Method of Intermediate Problems 

Let us suppose that .// is a separable Hilbert space 
with the inner product (u,v) and that A is a self- 
adjoint operator with domain (/ A in .//.' We will 
assume A is bounded below and that the first part 
of its spectrum is discrete. The eigenvalues and 
eigenvectors in this part of the spectrum will be 
denoted by X, and u v (v=l } 2, . . .) respectively; 
the first limit point of the spectrum will be called 
X ffi . Thus we have 



and 



Au v — \ v u p =0, (p=l,2, ...) 



Xl<X 2 <...<Xa 



(1) 



(2) 



Suppose further that A can be written as the sum 
of two operators, 



A=A°+A', 



(3) 



where A is a self-adjoint operator with a known 
spectral resolution E% and A* is symmetric and 
positive. The domains of A and A f are denoted by 
QP A and y A , respectively, so that £/> A =C/ A ,,n f j x . 
The positivity of A f is expressed by 



(A'u,u)>0 



(4) 
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for every u in the domain of A'. Finally we assume 
that the first part of the spectrum of A is discrete, 
and, listed in increasing order, the eigenvalues are 



X?<X2<... 



(5) 



The corresponding eigenfunctions are u\ 
Since A°<A, it follows 3 [4, p. 214] that 



< 



X?<X„ 



and 



X° m <\„ 



1 2 



(6) 



These inequalities go back to H. Weyl [8]. Thus 
the known eigenvalues of A give rough lower bounds 
for the initial eigenvalues of A. For this reason the 
operator A is called a base operator and the eigen- 
value problem for A is called a base problem. 

The procedure of intermediate problems links the 
base operator A to the given operator A by a sequence 
{A k } of intermediate operators. These operators sat- 
isfy the inequalities 



A°<A k <A k + l <A, (Jfc=l,2, . . .); 



(7) 



and consequently their eigenvalues give intermediate 
lower bounds to those of A. 

In order to define these intermediate operators, 
we introduce an auxiliary Hilbert space Jt '. M" ' is 
just the completion of Q) A > in the norm generated 
by the inner product [u,v] defined by 



[u,v]=(A'u,v). 



(8) 



For simplicity we will assume that A' is positive 
definite; the indefinite case requires only a slight 
extension of the theory. Let {p t } be a linearly inde- 
pendent set of vectors in Q) A >\ note that since A' 
is positive definite the sequence {A'pi} is also linearly 
independent. Define P h to be the projection in $?' 
on the span £/* of the first k elements of {pt}. For an 
arbitrary vector v in ,yf", P k v is given by 



PV=X) [v,Pi\1>tjPj, 



U=i 



(9) 



where the constants b tj are the elements of the 
matrix inverse to that with the elements {A!p u p d ). 
Since the subspaces Sf n are increasing with k; i.e., 



^c^c/, (A=l, 2, .. .), 



(10) 



the corresponding projections satisfy the inequali- 
ties, 



0<P k <P k + 1 <I, 



(11) 



where / is the identity operator in J^\ We now 
form the operator A'P k , which is well defined on 



3 A°<A means =^a 3 -& a .and (A°u,u)< (Au,u) for every u in -# a. 



Q) A . For an arbitrary vector w in Q A >, A'P k w is given 

by A'P k w=£ (v>*A'Pt)b tJ A'p„ (£=1,2, . . .). (12) 
u=i 

In the form given on the right in eq (12) A'P k has 
an obvious extension to all of J^f and there it is a 
symmetric completely continuous operator. Hence- 
forward we will regard A f P k as thus extended. 
It follows from eq (11) that in ^f the operators 
A r P k satisfy 'the inequalities, 

0<A'P k <A'P k+l <A f , (k=l, 2, . . .). (13) 

The intermediate operators A k are now defined b} T 

A*=A»+A'I*, (Jfe=l, 2, . . .); (H) 

k is called the order of the intermediate operator. 
Since the operators A'P k are symmetric and bounded, 
it follows [5, p. 301] that the operators A k are self- 
adjoint and have the same domain as A . And 
moreover, since A f P k is for each k a completely 
continuous operator, each of the operators A k has 
exactlv the same limit points in its spectrum as does 
A [5,"p. 367]. From eq (13) it follows directly that 

A°<A k <A k + l <A, (k=l, 2, . . .). 

These inequalities are just those stated in (7), and 
they show that the operators A k defined by eq (14) 
are in fact intermediate operators. This means that 
in the initial discrete part of the spectra the eigen- 
values satisfy the inequalities 



\° y <\ k <\ k+1 <\ v , (*=1,2, 



>), (15) 



where the superscripts on the eigenvalues indicate 
the order of the intermediate operator to which 
they belong. 

We now turn to the problem of determining the 
eigenvalues and eigenvectors of the intermediate 
operators A k . We will confine our attention to 
those eigenvalues which lie below the first limit point 
X°co common to all operators A k . 

The eigenvalue equation for A k , 



A k u-\u=0, 



(16) 



can be written in a more useful form by using the 
defining eq (14) for A k . This gives 



A°u-\u = -A'P k u, 
and by using eq (12) we obtain 



A°u— \u= 



:- 2 faA'pJbtjA'pj. 

i, 3= 1 



(17) 



(18) 



For convenience we designate the quantities y^.i 
(UjA'p^bij by a.j } so that the equations to be satis- 
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ficcl become 



and 



A°u— \u= — ^ajA'pj, 



(19) 



T J a J { Pj ^p l )={u^p l ), (Z=l,2, . . , k). (20) 



If u is to be an eigenvector corresponding to a 
value X which is not in the spectrum of A°, then u 
may be written 

u=-^2a J IhA'p j , (21) 

i=i 

where R y is the resolvent operator for -4°. 
The eqs (20) become 



^ajipj+Rui'PjA'pd^O, (Z=l,2, . . ., k). (22) 
j=i 

All of the a's cannot vanish if u is to be Qontrivial; 
therefore X must be such that the determinental 
equation 



\(p,+RxA%,A' Pi )\=0 



(23) 



is satisfied. 

Equations (22) have solutions for each value X that 
is not in the spectrum of A and that satisfies eq (23). 
The number of solutions corresponding to each such 
root is just the nullity of the coefficient matrix and, 
according to eq (21), these solutions give the eigen- 
vectors of A k associated with that root. These 
eigenvectors are clearly linearly independent since 
Rx is one-to-one on .//. 

For each value X that is an eigenvalue of the 
operator A°, say X°, the procedure is somewhat 
different. Let . //\ v be the characteristic subspace 
of A° at the eigenvalue X?. Suppose that it is 
spanned by the eigenvectors {v h v 2} . . . , Vi}. It is 
clear that if u is an eigenvector of A k corresponding 
to XJ, then in addition to the eq (19) and (20) the 
further restrictions, 

ZW^A-*) = 0, (i=k+l,k+2,...,k+l) (24) 
i-i 

must be satisfied. If we assume that eq (24) are 
satisfied, u is given by 

*=-2 aflyiA'Pi- 2 a^-», (25) 

;=i j=k+i 

where by R[o we mean the restriction of the resolvent 
operator for A to the part of our Hilbert space 
orthogonal to ^o. The eqs (20) become 

ft fc+7 

Jla j (p j + R^A'p j ,A'p i )+ J^ afa^A'p^O, 

j=l j=k-\-l 

(i=l,2, ...,fc). (26) 



These, together with eq (24), determine the constants 
oil. The necessary and sufficient condition for the 
existence of nontrivial solutions is the vanishing of 
the determinant of the coefficients; i.e., 

/(pj+R^oA'p^A'pi) ! (vj. k , A'pi)\ 
det I — 1=0. (27) 

\ (A'p^Vi-z) / 

Let us assume that (27) holds; that is, that X? is an 
eigenvalue of A k as well as of A . The eigenvectors 
of A k corresponding to XJ are then determined in the 
following way: The eqs (24) and (26) have linearly 
independent solutions {a\, cti, . . . , a° k+l }, (<r= 
1,2, . . . ,n), where n is the nullity of the coefficient 
matrix displayed in (27). By (25) these give rise 
to n vectors u a that are eigenvectors of A k for the 
value XJ. These are also linearly independent, for if 
C a are constants such at 2, = 1(7^=0, then by the 
expression (25) for u a we must have 

\<r = l j=\ / <7 = 1 J-k+1 

(28) 

But since R' x a is a one-to-one operator on the orthog- 
onal complement of , ^n, eq (28) implies the separate 
statements, 

i;(sa«jW,=o, g (saaj)»,_*=o. 

i=l \er=l / /=ft+i \«r=] / 

(29) 

The linear independence of t ho vectors A'pi and of 
the vectors Vj imply that 



s^ 



=o f 



(j=l,2,« ••,*+«). 



But since the vectors {a*, a° 2 , . . . , a£+z} arc 
linearly independent, each of the constants C must 
vanish, and thus the vectors u a arc linearly inde- 
pendent. 

As we observed previously, the part of the 
spectrum of A k that lies below X° is discrete. The 
eigenvalues there are completely determined from 
eqs (23) and (27), and the eigenvectors are found 
from corresponding solutions of the linear equations. 
These eigenvalues give intermediate lower bounds 
to those in the first part of the spectrum of A. 

For the sake of completeness and later use, we 
sketch a proof of convergence of the eigenvalues of 
the intermediate problems. The conditions that we 
use are those stated by Aronszajn [1]. We assume 
that A and A have completely continuous inverses, 
that A' is bounded relative to A , and that the set 
{pi} is complete in jf '. For convenience we will 
assume that A is positive definite; the modifications 
for A bounded below are straightforward. 

The purpose of the proof is to show that the inverse 
operators converge uniformly to A' 1 . From this 
the convergence of the eigenvalues follows immedi- 
ately [5, p. 372]. 



5SG066— 61- 
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Since the inverse operators (A*)' 1 form a decreas- 
ing sequence of symmetric operators bounded below 
by A~ l , the sequence converges strongly to some 
symmetric operator [5, p. 263]. The strong con- 
vergence to A -1 is guaranteed by the inequalities 
that hold for any t in J% 

\\[(A«)-i-A^]tf 

^IK^ )" 1 -^- 1 !! 3 !!^!!^^)- 1 -^- 1 ]^^), (30) 

and 

{[{A k )- l -A- l ]t,ty<[{I-P k )A-H,{I-P k )A-H] 

.(A'(A k )-H,(A k )-H). (31). 

The first of these inequalities follows by using the 
Schwarz inequality in the form, generated by (^4 fc ) _1 — 
A' 1 and then using the fact that (A*) _1 < (A°)-\ 
The second is obtained by writing (A k )~ 1 —A~ 1 = 
{A k )- l A'(I-P k )A-\ and then using the Schwarz 
inequality in the form generated by A' . The first 
factor on the right in (31) vanishes with increasing 
k since the set {p t } is complete in jf f . The second 
remains bounded since A' is bounded relative to A . 
The strong convergence of (^L*) -1 to A' 1 is thus in- 
sured. However, since each operator (A k )~ l as well 
as A~ l is completely continuous and symmetric, the 
mono tonic strong convergence implies uniform, con- 
vergency by an easily proved analog of Dini's 
theorem. 

Remark: We observe that if the span of the 
vectors {#*}*=! contains some eigenvectors [u v , u Vt2 
. . . , u v ] of A, (j<k), then A k =A on the span of 
these u's so that X„ , X v . . . , X„. are eigenvalues of 
A k as well as A, and they have the corresponding 
eigenvectors u vv u V2 , . . . , u v ., respectively. In- 
vestigations based on this property lead to state- 
ments of the following type: 

Theorem: Let A be an operator with a purely 
point spectrum, and let X„ be an eigenvalue of A that 
is less than Xj?, then there exists an intermediate 
problem, (of sufficiently high order) which has 
XJ=X M and u k =Up for ju=l,2, . . . , v. 

3. Truncation of the Base Operator 

The principal difficulties in the computation of 
approximating eigenvalues and eigenvectors by the 
method just described arise from the fact that the 
resolvent operator R\ for the base problem rarely 
is known in closed form. In most cases it can be 
expressed only by use of infinite sums or integrals; 
i.e., in our cases by 



n\W—2-j "To — x 



r* oo 

Jx . 



dE°w i 

}JL — X 



We overcome this difficulty by changing the base 
problem so that it has a simple resolvent and still 
provides lower bounds by the method of inter- 
mediate problems. These bounds are computed 
from equations identical with those given earlier ex- 



cept that the resolvent is expressed by a finite sum 
with a remainder in closed form. 

Let us assume that A has infinitely many eigen- 
values less than >«, (the modifications for the cases 
in which there are only a finite number will be im- 
mediate). Define the truncated operators A i>0 by the 
equations 

A^=A Q E^+\hi(I-Elo), (t=l,2,...), (32) 

and denote the corresponding resolutions of the 
identity by E{'°. It is clear that each operator A i>(> 
is bounded and symmetric, has the eigenvalues and 
eigenvectors X?, u° for p=l,2, . . ., i, and also has 
\ i + ^ as an eigenvalue of infinite multiplicity. The 
resolutions of the identity may be written down im- 
mediately. They are 



K - U x? +1 <x/ ; 



(i=l,2,...)- (33) 



An operator A*' is called a truncation of A of order i. 
These operators satisfy the fundamental inequalities 

A^KA^^KA^ (i=l, 2, . . .), (34) 

since from the spectral theorem it follows that 

([^+ 1 '°-^°]^,u)=(XS +2 -X? + i) r _ d(E x u,u)>0. 

for each u in J% and that 

(Ui°-A''°Ktt)= f (X-X? +1 )d'(#x^)>0 



for each u in Q} A *. 

The new intermediate operators A uk based on 
A i>0 are defined by 

A i ' k =A i '°+A'P k , (i,k=l,2, . . .). (36) 

These are bounded, symmetric, and are monotoni- 
cally increasing in i and k; i.e., 



A i > k <A i+1 > k <A k <A, (i,k=l,2, . . .), (37) 



and 

A^KA'^KAW+A'KA, (i,*=l,2, . . .). 

(38) 

The montonicity in the first index follows from (34) 
and in the second index from (13). The inequalities 
(37) and (38) imply the parallel inequalities for the 
corresponding eigenvalues, 



y<\l +i ' k <\ k v <\ v | 



x^<xj»* +1 <x„ 



(,= 1,2,...), (39) 
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where by Xj ■* we mean the i>th eigenvalue of A*'*. 

Since A i,k and ^4 z0 differ by a completely con- 
tinuous operator, the limit point X? +1 persists and 
is the only limit point of the spectrum of the inter- 
mediate operator A Uk as the order k of the operator 
is increased [5, p. 367]. Further, the inequalities 
(39) show that for any k there are at most { eigen- 
values of A Uk less than X? +1 . 

The spectrum of each operator A i<k can be deter- 
mined easily by the usual procedure for intermediate 



problems, 
stated by 



The eigenvalue equation for A i,k is 



A ik u—Xu=0, 



or equivalently (in the notation of the preceding 
section) 



A i -°u-\y = — '22a j A'pj, 



If u is an eigenvector corresponding to a value X 
that is not in the spectrum of A 1 ' , then u may be 
written 






(40) 



where R£ is the resolvent operator of A i>0 . We 
note here that for any vector v in 2% R{v is given 
by the closed expression, 






The eqs (22) become 



sj^^S^^^} 



(41) 



^^(p j +R{A%,A'p l ) = J (J=l,2,. •-,*). 

(42) 

As a consequence, the necessary and sufficient condi- 
tion for such a X to be an eigenvalue is that it be a 
root of the equation 



\( Pi +RlA' Pj ,A' Pl )\ = 0. 



(43) 



The multiplicity of each such eigenvalue is the nullity 
L of the coefficient matrix of eq (42). The eigen- 
vectors of A Uk are determined from the solutions 
of eq (42) by eq (40). 

Now suppose that \° v is an eigenvalue of A f '° that 
is less than X? +1 and that the characteristic sub- 
space M\ v of A f '° corresponding to XJ is spanned by 
[V], v 2 , . . , Vi}. If u is to be an eigenvector of 



A i>k corresponding to Xj, then u is given by the 
equation 



u~- 



-2>, h% A'pj- §] ajVj_ k , (44) 

7=1 ;'=A-+i 



where the constants aj must satisfy the equations 



^afa+B$A'p^'p r ) 

3=1 

+ ?& afa-*A%)=0, (r=l,2,...,i) (45) 
and 



Sa,(^^V-,) = 0, (r=k+l,k+2,...,k+l). (46) 
;=i 

The linear independence of the eigenvectors arising 
from Linearly independent sets of a's may be demon- 
strated as in the preceding section. 

The multiplicity of X 1 ,! as an eigenvalue of A ulc is 
thus just the nullity of the matrix 



fipj+R&A'p^A'pr) ! (Vj- k ,A'p r )\ 



(47) 



(A'pj, v r _ k ) 







We have found all of the eigenvalues of A Lk that 
are not equal to X? +1 . The complete continuity of 
A'P k implies that X? +1 , the only limit point of the 
spectrum of A it0 , remains as the only limit point of 
the spectrum of A uk . Thus the total spectrum of 
A uk is determined. The spectral subspace of A itk 
associated with X ( - +1 is just the orthogonal comple- 
ment of the span of all the eigenvectors already 
determined for the values of X not equal to X? +1 . 

Thus, as stated at the beginning of this section, 
the method of truncation justifies the replacement 
of the base problem resolvent by a suitable approx- 
imation in closed form. The necessary computa- 
tions can be carried out by routine techniques. 

The convergence of the eigenvalues of the trun- 
cated intermediate operator A i>k to those of A k can 
be shown under the assumption that A and A have 
completely continuous inverses. We demonstrate 
the uniform convergence of (A*-*) -1 to (A* )' 1 for 
each fixed k. For convenience we will assume that 
A° is positive definite. First we note that the 
operators (A*' )' 1 converge in norm to (A )' 1 as i 
becomes infinite. This follows directly from the 
spectral theorem, for we have 



(A°y 



-r 



\-UE$, 



and 



(A'-°)- ,= 



/ 






A i+iJxL-0 



dEl, 
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so that 



Consequently we have 

\\(A")- 1 -(A i -oy i \\<4-' (48) 



•(^•°)- 1 ||< I f- 



Since \° iVi becomes infinite with i, the norm con- 
vergence of (^4'°) _1 is clear. This implies that 
CA M ) -1 converges in norm to (A*)' 1 . The proof 
goes as follows. We start from the simple identity 

A i ' k -A k =A i '°-A°. 

From this it follows immediately that 
4«.^(4*)-i-( i l<.*)-iJ i l*= - 4<.o[(^o)-i-( 2 li.o)-i] i io j 

and thus 

(^)" 1 -(^'' fc )" 1 = (^ , ")" 1 ^'° 

[(^L )- 1 -^ )- 1 ]^! ^)- 1 . (49) 

By writing (A*'*)" 1 ^ as I—(A i ' k )- 1 A'P k it is evi- 
dent that (.A* ,fc )"M.*' is an operator that can be 

bounded independently of i since IK^l 1 *) -1 !! <ro* 

1 
Since A®(A k )~ l is also bounded for each fixed k, the 

convergence in norm of (A itk )~ l to (A k )~ l follows 

directly from (48) and (49). 

The uniform convergence of (^4 f,A ") _1 to (A k )~ l 

and that of (A k )~ l to A~ l demonstrated in section 2 

are sufficient to insure that the eigenvalues of A itk 

converge to those of A as i and k become large. 

4. Truncation by Choice of Elements 

In this section we discuss an extension of another 
procedure [2] for truncating the resolvent R\ that 
appears in section 2. This can be accomplished by 
choosing the elements p t in suitable ways instead of 
changing the base problem operator as has been 
done in previous sections. 

We assume that the vectors {p u p 2 , . . ., p k ] which 
determine A k have been chosen so that 



N 



A f p t =J2 ft,< (i=l, 2, . . .,fc). 



(50) 



That is, it is assumed that each of quantities 
A'pi(i=l, 2, . . ., k) can be expressed in terms of a 
finite number N of eigenvectors of A . The number 
N will of course depend on k. The quantities 
(R\A'pi,A'pj) which appear in section 2 will have 
the simple expressions 



(BiA' Pi ,A%)=i; $$2%, 



and thus are effectively truncated by the choice of 
elements (50). The computations for the eigen- 
values and eigenvectors ma}' be carried out directly 
from the equations of section 2. 

In the present case, however, it is simpler to 
proceed directly. We have already shown that the 
expression A k u can be written 



A k u=A°u+ 2 (%A'pi) bijA'pj, 



where, as before, the quantities 6^ are the elements 
of the matrix inverse to that having (A'pi,p } ) as 
elements. Using (50) this expression becomes 

A k u=A°u+^b S M)foAA*< (52) 

i,j=l n,v=l 

We can now obtain all of the eigenvalues and eigen- 
vectors of A k by inspection. First consider those 
vectors in ^ A o that are orthogonal to every eigen- 
vector of A . If u is such a vector, then (52) shows 
that A k u—\u=A°u—\u. Clearly no such a vector 
can be an eigenvector of A k . Further, if u is an 
eigenvector ul of A with an eigenvalue A° and is 
not used in the sums in (50), then by (52) we have 
A k ui—\° ff u a =A u° a — \ ff u a =0. Hence each such 
vector is an eigenvector of A k with the same eigen- 
value as for A° ff . 

It remains to consider those vectors u which 
are linear combinations of the vectors u° v (v=l, 
2, . . ., N) used in the sums in (50). There are 
exactly N linearly independent vectors of this kind 
that are eigenvectors of A k . These are determined 
as follows. We write u as 



u^ 



N 



(51) 



and insert this in the eigenvalue equation 
A k u— \u=0, where A k is given by (52). The 
linear independence of the vectors u° v leads to the 
equivalent algebraic system, 



St/s M<Am+(x!-a)0 =0, 

0i=l,2, . . ., N). 



(53) 



Consequently the determination of the intermediate 
eigenvalues and eigenvectors depends only on the 
solution of a matrix eigenvalue problem. 

This procedure can be especially useful when the 
operator A f is a multiplicative operator, as for 
example in the spheriodal wave equation. 4 In the 
prolate case the eigenvalue equation is stated by 

Au— \u=0, 



4 The numerical applications in this and other cases are presently being carried 
out by W. Bbrsch-Supan. 



no 



in which A is the essentially self-adjoint operator 
given by 



A d 

Au= — r- 
dx 



[ 



(1- 



J dxyi 



///- 



;U-\-C 2 X 2 U 



on the set of all functions continuous on I — 1, +1] 
that are twice continuously different iable on 
(—1, +1) and for which An belongs to J5^ 2 (— 1, +1). 
Here m is a non-negative integer and c is a real 
constant. A base operator is obtained by omitting 
the positive term A' = c l x 2 . The solutions of the base 
problem are associated Legendre functions. The 
choice pi=u° i =F h Z l+ i-i and the recursion relation, 
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2n+l 



allow one to express A'pi as a linear combination of 
u°i- 2 , w$-i, u° if u° i+ i, and u? +2 . Other examples of 
this sort can be found easily. 
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